Exact Solution of Strongly Interacting Quasi-One-Dimensional Spinor Bose Gases 
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We present an exact analytical solution of the fundamental system of quasi-one-dimensional spin-1 bosons 
with infinite 5-repulsion. The eigenfunctions are constructed from the wave functions of non-interacting spinless 
fermions, based on Girardeau's Fermi-Bose mapping, and from the wave functions of distinguishable spins. We 
show that the spinor bosons behave like a compound of non-interacting spinless fermions and non-interacting 
distinguishable spins. This duality is especially reflected in the spin densities and the energy spectrum. We find 
that the momentum distribution of the eigenstates depends on the symmetry of the spin function. Furthermore, 
we discuss the splitting of the ground state multiplet in the regime of large but finite repulsion. 
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Introduction. — Strong correlations are the basis of some of 
the most fascinating and important phenomena in many-body 
quantum systems. One intriguing only recently experimen- 
tally realized example HHH is the Tonks-Girardeau gas J^, 
in which a one-dimensional (ID) system of strongly interact- 
ing bosons acquires fermionic properties. Already in 1960 
Girardeau found an elegant solution of these systems which 
maps the wave functions of non-interacting spinless fermions 
to that of spinless hard-core bosons. However, a thorough 
description of quasi- ID bosons needs to include the spin de- 
grees of freedom. Surprisingly, thus far no Fermi-Bose map- 
ping for these important systems exists. In current investiga- 
tions Girardeau's Fermi-Bose mappinghas been extended to 
fermionic Tonks-Girardeau gases H-SSHland ver Y recently 
also to mixtures [8] and two-level atoms |9J]. Ferromagnetic 
behavior in a strongly interacting two-component Bose gas 
has been found in the thermodynamic limit 11 1 Oh . 

In this Letter we present for the first time an analytically 
exact solution of quasi- ID bosons with infinite J-repulsion, 
which includes the spin degrees of freedom. This opens up 
new possibilities for a precise manipulation and description 
of these fundamental systems. We show that the strongly 
interacting spinor bosons behave like a compound of non- 
interacting spinless fermions and non-interacting distinguish- 
able spins which is especially reflected in the spin densities 
and the energy spectrum. These properties are very differ- 
ent from those of weakly interacting spinor Bose gases 11 ill . 
Another key result is that in contrast to single-component sys- 
tems the momentum distribution of the spinor bosons depends 
on the symmetry of the spin function. 

Model. — The many -body Hamiltonian describing A quasi- 
1D spin-1 bosons at zero magnetic field is given by the 
3 N dimensional matrix 111 111 



H = 



i = \ L 1 

+ ^2s(zi - z-j) 



H — raw 2 ?,- 



U t spin + U 2 fi-f j 



(1) 



where m is the mass of the bosons, uj is the axial trap fre- 
quency, 11 S pi n is the identity matrix in spin space, Uo and U~2 
are the coupling constants of the spin-independent and spin- 
dependent interactions, and /j denote the dimensionless spin- 
1 matrices of boson i. We consider a highly elongated trap 
where the transverse trap frequency ui± is much larger than the 
axial trap frequency (u>± 3> w) so that the transverse motion 
is restricted to zero-point oscillations. In this limit the effec- 
tive coupling constants are given by Ui = 2huj±Cij (i = 0, 2) 
with 7 = 1/(1 - 1.4600^/210 |Q, c = (ao + 2a 2 )/3, 
and C2 = (ct2 — ao)/3 111 ill , where ao and ai are the s-wave 
scattering lengths and l± = y/h/(mu)±). We neglect p- 
wave scattering [6]. Thus, Uq and U~2 are proportional to the 
transverse trap frequency u±_ which has been used in experi- 
ments [2] to reach the strong interaction limit. The Hamilto- 
nian (UJ conserves parity n, total magnetization F z , and to- 
tal spin F 2 (F = fi)- In the limit of infinite repulsion 
(Uo = oo) we use a Fermi-Bose map to construct the eigen- 
functions of (Q]i, whereas for large but finite repulsion we nu- 
merically diagonalize ([TJ via a finite basis set approach lfl3ll . 

Exact eigenstates in the limit of infinite repulsion. — We 
first construct the exact eigenfunctions of the Hamiltonian (UJ 
in the limit of infinite repulsion (Uq = oo). In this limit one 
can neglect the spin-dependent interaction and we set U-2 = 0. 
For simplicity we restrict our discussion to the ground 
states of spin-1 bosons. Because of the infinite <5-repulsion 
the configuration space M. N decomposes into A! sectors 
Cir = {(zi, ...,z N )£ R N , z^x) < • • • < z^fff)}, where ir 
is an arbitrary permutation. At the boundaries of these sec- 
tors the wave function has to be zero [ip(zi, . . . , Zn) = if 
Zi = Zj] and within these sectors the wave function has to 
obey the Schrodinger equation of non-interacting particles. 

Let us disregard the spin of the wave function in a first step. 
The Fermi ground state tp^ of the problem above is given by 
the ground state Slater determinant of non-interacting spin- 
less fermions. The ground state of spinless bosons ipE K 
given by the absolute value of the Fermi ground state [3]: 
0Q 3 = \ipo\. There is no symmetry restriction to the wave 
function of distinguishable particles. Thus, any restriction of 
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the Tonks-Girardeau ground state ip^ to an arbitrary sector CV 
is a non-symmetric ground state of JV spinless distinguishable 
particles which we denote by |tt). Restrictions to different 
sectors are orthogonal. Thus, the ground state of JV spinless 
distinguishable particles is AH times degenerate. An orthonor- 
mal basis of the space of ground states is therefore given by 

I \ J ko] if < ■ ■ ■ < Z <N) ,~ 

(zx, . . . , z N \n) = < , (2) 
10 otherwise. 

Now we take the spin into account. The product of an ar- 
bitrary orbital function \ir) with an arbitrary TV-particle spin 
function \x) is an eigenfunction of Hamiltonian CO, since H is 
diagonal in spin space (we have neglected U^). Thus, the wave 
function |7r)|x) describes distinguishable hard-core particles 
with spin. In order to describe bosons one has to symmetrize 
this non-symmetric wave function. One can directly construct 
all bosonic solutions from the spin functions by means of the 
unitary map 

W\x) = VWPs(\id)\x)), (3) 

where Ps = ^7 Su-eSiv ^( 7r ) denotes the symmetric projec- 
tion operator and id is the identical permutation. The action 
of the permutation operators U(ir) onto the spinor wave func- 
tions is given by 

U(TT)\TT')\fi,...,f N ) = |7r7r')|/ n -i(i),---,/ 7r -i(iv))> 

where tttt' is the composition of the permutations it and it', 
and fi = ±1, are the spin indices of the ith particle. Equa- 
tion ([3J, as the main result of our Letter, provides an analyt- 
ical and general extension of Girardeau's Fermi-Bose map- 
ping to bosons with spin degrees of freedom. The map W 
is isometric, and therefore, an orthonormal set of spin func- 
tions is directly mapped onto an orthonormal set of bosonic 
ground states. It is also bijective, and thus the degeneracy of 
the ground state equals the dimension of the JV-particle spin 
space. Further, we can directly construct bosonic F z and F 2 
eigenfunctions VF|x)> since F z and F 2 commute with Ps- 

The above construction scheme can be extended to the ex- 
cited states and to bosons with arbitrary spin. To obtain the ex- 
cited states one has to replace | | in Eq. (O by Aijjf, where 
ijjf is the ith eigenfunction of the non-interacting fermionic 
problem and A is the "unit antisymmetric function" yfl. It 
follows immediately that the energy eigenvalues of the spin-1 
bosons at zero magnetic field agree with the energies of spin- 
less non-interacting fermions; i.e., the ground state energy is 
given by E g = N 2 /2 frio and the level spacing is AE = 1 tuo. 
However, the degeneracy of each energy level is 3 N times 
larger than in the spinless case. A homogeneous magnetic 
field lifts the degeneracy of the levels and shifts the energy of 
each F z eigenstate W\fi, . . . , /at) according to the Zeeman 
energy of the spin function . . . , /jy). 

Large finite repulsion. — In the following we analyze the 
structure of the ground state multiplet in the regime of large 
but finite repulsion. The results are based on a numerical diag- 
onalization of ((T). In the limit of infinite repulsion, Uq = 00, 
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FIG. 1: (color online) Cut through the non-zero spin-components 
of a bosonic 3-particle state (see text). The second and the third 
coordinate are fixed at Z2 = —0.8 1 and 23 = 0.8 /. Shown is the 
exact wave function in the limit of infinite 5-repulsion (red solid line) 
and the solution of a numerical diagonalizion of (T) in the limit of 
large but finite repulsion (blue dashed line). 

the ground state is 3^ times degenerate. For realistic sys- 
tems the spin-independent interaction Uq is large and finite 
and the spin-dependent interaction U2 is much smaller so that 
the ground state energy level is quasidegenerate. Let us first 
consider the case that U% is zero. Then, states with completely 
symmetric spin functions have lowest energy lfl4ll whereas 
states with most antisymmetric spin functions have highest 
energy. The energy gap AEo, which arises from the sym- 
metry of the spin function (which is closely related to the 
symmetry of the corresponding orbital functions), is propor- 
tional to Hlu/Uo \Ui = Ui/(hwl) are dimensionless interac- 
tion strengths and I = \J h/(muj ] . A similar energy structure 
has been discussed in Ref. [8] for a two-component system. 
Let us now switch on 1/2- Then, all states become F 2 eigen- 
states and the states with completely symmetric spin functions 
have total spin F = N, N - 2, N - 4, ... lUj. For fer- 
romagnetic coupling, U2 < 0, states with F — N have low- 
est energy whereas for antiferromagnetic coupling, U2 > 0, 
states with F = 1 (if JV is odd) or F — (if JV is even) be- 
come ground states. The energy gap AE 2 , which arises from 
the spin-dependent interaction, is proportional to hui U 2 /Uo, 
since the local correlation function g 2 is proportional to X/Uq 
in the limit of strong repulsion [ 15] and AE 2 = ^52- 

For a sufficiently large repulsion, the quasidegenerate 
ground states are well approximated by the states (O. As an 
example, Fig. Q] shows the 3-boson state \ip) = W\x) with 
| X ) = 1/2(|0,1,-1) + |0,-1,1)-|1,-1,0>- I -1,1,0)) 
(red solid line), which is a F z , F 2 , and parity eigenstate 
(F z = 0, F = 1, II = — 1), compared to the corresponding 
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solution of a numerical diagonalization of (fl} with parameters 
U /l = 20 hto and U 2 = -J7 /2000 (blue dashed line). Both 
solutions are in excellent agreement, and thus the limit of infi- 
nite repulsion is practically reached at Uq/1 = 20 fko lfl3l[l6ll . 
Note that at finite Uo the state shown in Fig. [TJ is an excited 
state within the ground state multiplet, since |x) is not com- 
pletely symmetric and F < 3. 

Spin densities of the ground states. — Apart from the energy 
spectrum, the dual nature of the hard-core bosons with spin 
can best be seen in the densities. While the total density is 
independent of the spin degrees of freedom and equal to that 
of non-interacting fermions, the spin density 



Pfz 0) 



E 



J dzi...dz N ^2 S(z-Zi)5f zfl 



x\^fu...J N (zi,---,z N )\ 2 



resembles the one of a chain of localized, distinguishable 
spins. In the limit of infinite 5-repulsion the wave function 
is given by \ip) = W\x) and one finds 



(4) 



with the probability pi (/ z ) that the zth spin is equal to f z 

/'-(/,)= >; \(A,---jN\x)\ 2 Sf, fi , 



E 



and the probability density pW(z) to find the ith particle of 
the system, restricted to the standard sector Ca, at point z, 

P {l) (z) = N\ f dz 1 ---dz N \^(z 1 ,...,z N )\ 2 

J 21 <-<Zjv 

xS(z ~ Zi). 



An explicit calculation yields the following formula 



(0 = d (*" {-l)»-*(N-k-l)\ 
dz {i - 1)\(N - k - i)M 



x^det[B(z)-Al]| A =o J, 

where the N x A-matrix B(z) has entries /3ij(z) — 
J_ dx4>i(x)<pj(x) with the single-particle eigenfunctions of 
the spinless problem </>j. 

Spin densities of selected ground states of 8 spin-1 bosons 
are shown in Fig. [2] Since the total spin of each atom is 
one, the spin density has 3 components which correspond 
to f z = ±1, 0, and which are drawn as a blue dashed (pi), 
red solid (po), and a green dotted line (p-i). By interpret- 
ing the square root of the densities pW (gray dash-dotted line 
in Figs. [2^ and b) as the wave packet of an imaginary par- 
ticle, : — y / pM", we obtain the intuitive picture of dis- 
tinguishable localized particles which are arranged in a row 
along the z-axis, one after the other |0, fl7tl . Each particle 
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FIG. 2: (color online) Spin densities of 8 spin-1 bosons in different 
ground states (see text). Shown are the densities p' ! ' (gray dash- 
dotted line, see text), and the components po (red solid line), pi (blue 
dashed line) and p_i (green dotted line) of the spin density. The spin 
densities resemble the one of a chain of localized spins. 



has a certain spin orientation fa according to the spin state 
Ix) = l/i, hi • • ■)■ The spin density of the bosonic system 
now corresponds to the spin density of this row of distinguish- 
able particles with the spin orientations given by |%). 

Fig. [2h shows the spin density of the spin-polarized state 
= 1) ■ ■ ■)■ 1° tne spin-polarized case all the 

probabilities p%{l) = 1 and Eq. (0]) reduces to: p\ = 
Si = PFermi- Fig. [2t> shows the spin density of the 
ground state W\l, 1, 1, 0, 0, 0, 0, 1). Similarly one has to add 
the particle densities p^ — p^ and p^ to the component 
pi, and the particle densities p^ — p' 7 ' to the component 
po of the spin density. Finally, Figs. [2t and d show the spin 
densities of the ground state W\l, —1, 1, 0, 0, — 1, — 1, 1) and 
the superposition state W-^ (|0,0, . . .) + |1, -1, 1,-1, • ■ ■)) 
respectively. 

Momentum distributions of the ground states. — One of the 
most important experimentally accessible quantities is the mo- 
mentum distribution of the spinor bosons, given by 

p(p) = E J d Pi--- d PN E 5{ <p ~ p ^ 

x IV'/i,...,/ n (pi,---,Pjv)| 2 - 

Fig. [3] shows selected momentum distributions of 5 spinor 
bosons in their degenerate ground states, obtained from a 
numerical diagonalization of dTJ. For comparison we have 
also plotted the momentum distribution of 5 non-interacting 
fermions (gray dashed line). It turns out that the shape of 
the momentum distribution depends on the symmetry of the 
spin function so that different ground states can have com- 
pletely different momentum distributions. States with a com- 
pletely symmetric spin function have a completely symmet- 
ric orbital function, VF|xs) = IV'ol |x«)> which is given by 
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the symmetry of the spin functions. Some ground states have a 
momentum distribution which is much broader and flatter than 
for a spinless Tonks-Girardeau gas. For large but finite repul- 
sion we have discussed the level structure of the ground state 
multiplet. Due to the similarity with distinguishable spins the 
spinor Tonks-Girardeau gas might be a promising candidate 
for novel quantum computation schemes. 
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FIG. 3: (color online) Momentum distributions of 5 spin-1 bosons in 
different ground states. The gray dashed line shows for comparison 
the momentum distribution of 5 non-interacting fermions. The shape 
of the momentum distribution depends on the symmetry of the spin 
function. 



the usual Tonks-Girardeau wave function. The momentum 
distribution of these states is equal to the one of a spinless 
Tonks-Girardeau gas (blue dash-dotted line) which exhibits 
a pronounced zero-momentum peak and long-range, high- 
momentum tails 11131 llal . The other extreme case is given 
by a flat and broad momentum distribution which resembles 
the fermionic one (red solid line). In the case of 2 and 3 spin- 1 
bosons some spin functions |% a ) can be completely antisym- 
metric and thus the corresponding ground state is given by 
W|Xa) = ipo\Xa) so that its momentum distribution is equal 
to the one of spinless fermions. One cannot construct com- 
pletely antisymmetric spin functions with more than 3 spin-1 
particles. However, it is possible to construct non-symmetric 
spin functions which are "almost antisymmetric" (see Young's 
Tableaux 01911 ) resulting in momentum distributions which are 
almost fermionic. We believe, due to group theoretical argu- 
ments [19], that this broadening and flattening of some mo- 
mentum distributions saturates for large N. 

Preparation and manipulation of the states. — Spin polar- 
ized Tonks-Girardeau gases have been prepared with ultracold 
87 Rb atoms in optical lattices JH 01 • To prepare states with 
non-symmetric spin functions we suggest to apply a radio fre- 
quency together with a magnetic field gradient so that the spin 
of each atom is rotated into a different final state depending 
on its position. 

Summary. — We have constructed the exact eigenstates of 
the fundamental system of quasi- ID spin-1 bosons with in- 
finite (5-repulsion by means of a Fermi-Bose map. The con- 
struction scheme and the formula for the spin densities are 
valid independent of the particle number, the spin of the 
bosons, and the (spin-independent) trapping potential. From 
the exact solutions we have determined the energy spectrum 
and the spin densities of the ground states. The momentum 
distribution of the eigenstates turned out to be dependent on 



* Electronic address: fdeuretz@physnet.uni-hamburg.de 
t Electronic address: klaus.fredenhagen@desy.de 

[1] B. Paredes etal, Nature (London) 429, 277 (2004). 

[2] T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305, 1125 
(2004). 

[3] M. D. Girardeau, J. Math. Phys. 1, 516 (1960); Phys. Rev. 139, 

B500 (1965), Sees. 2, 3, and 6. 
[4] T.Cheon and T.Shigehara, Phys. Lett. A243, 111 (1998); Phys. 

Rev. Lett. 82, 2536 (1999). 
[5] M. D. Girardeau and M. Olshanii, Phys. Rev. A 70, 023608 

(2004). 

[6] B. E. Granger and D. Blume, Phys. Rev. Lett. 92, 133202 
(2004) 

[7] S. A. Bender, K. D. Erker, and B. E. Granger, Phys. Rev. Lett. 

95, 230404 (2005). 
[8] M. D. Girardeau and A. Minguzzi, Phys. Rev. Lett. 99, 230402 

(2007). 

[9] S. V. Mousavi, A. del Campo, I. Lizuain, and J. G. Muga, Phys. 

Rev. A 76, 033607 (2007); M. D. Girardeau, arXiv:0707.1884. 
[10] X.-W. Guan, M. T. Batchelor, and M. Takahashi, Phys. Rev. A 

76, 043617 (2007). 
[11] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998); T. Ohmi and 

K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998); C. K. Law, 

H. Pu, and N. P. Bigelow, Phys. Rev. Lett. 81, 5257 

(1998); J. Stenger et at, Nature (London) 396, 345 (1998); 

M. R. Matthews et al, Phys. Rev. Lett. 81, 243, (1998). 
[12] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[13] F. Deuretzbacher, K. Bongs, K. Sengstock, and D. Pfannkuche, 

Phys. Rev. A 75, 013614 (2007). 
[14] E. Eisenberg and E. H. Lieb, Phys. Rev. Lett. 89, 220403 

(2002) . 

[15] D. M. Gangardt and G. V. Shlyapnikov, Phys. Rev. Lett. 90, 
010401 (2003). 

[16] Y. Hao, Y. Zhang, J. Q. Liang, and S. Chen, Phys. Rev. A 73, 

063617 (2006); S. Zollner, H.-D. Meyer, and P. Schmelcher, 

Phys. Rev. A 74, 053612 (2006). 
[17] I. Romanovsky, C. Yannouleas, and U. Landman, Phys. Rev. 

Lett. 93, 230405 (2004); O. E. Alon and L. S. Cederbaum, Phys. 

Rev. Lett. 95, 140402 (2005). 
[18] M. D. Girardeau, E. M. Wright, and J. M. Triscari, Phys. Rev. A 

63, 033601 (2001); T. Papenbrock, Phys. Rev. A 67, 041601(R) 

(2003) ; S. Zollner, H.-D. Meyer, and P. Schmelcher, Phys. Rev. 
A 74, 063611 (2006). 

[19] J. J. Sakurai, Modern Quantum Mechanics (Addison- Wesley, 
Reading, MA, 1994), Chap. 6.5. 



